Optimal Orientation of Anisotropic Solids 

A. N. Norris 

Mechanical and Aerospace Engineering, Rutgers University, 
Piscataway NJ 08854-8058, USA norris@rutgers.edu 

Abstract 

Results are presented for finding tiie optimal orientation of an anisotropic elastic 
material. The problem is formulated as minimizing the strain energy subject to rota- 
tion of the material axes, under a state of uniform stress. It is shown that a stationary 
value of the strain energy requires the stress and strain tensors to have a common 
set of principal axes. The new derivation of this well known coaxiality condition uses 
the 6-dimensional expression of the rotation tensor for the elastic moduli. Using this 
representation it is shown that the stationary condition is a minimum or a maximum 
if an explicit set of conditions is satisfied. Specific results are given for materials of 
cubic, transversely isotropic (TI) and tetragonal symmetries. In each case the exis- 
tence of a minimum or maximum depends on the sign of a single elastic constant. The 
stationary (minimum or maximum) value of energy can always be achieved for cubic 
materials. Typically, the optimal orientation of a solid with cubic material symmetry 
is not aligned with the symmetry directions. Expressions are given for the optimal 
orientation of TI and tetragonal materials, and are in agreement with results of Rovati 
and Taliercio I15| obtained by a different procedure. A new concept is introduced, the 
strain deviation angle, which defines the degree to which a state of stress or strain 
is not optimal. The strain deviation angle is zero for coaxial stress and strain. An 
approximate formula is given for the strain deviation angle which is valid for materials 
that are weakly anisotropic. 

1 Introduction 

The strain energy of a piece of homogeneous anisotropic elastic material depends on the 
orientation of the material relative to the directions of principal stress, although the ori- 
entation dependence vanishes trivially for isotropic solids. This property is therefore an 
inherently anisotropic feature of elasticity, and it raises the question of how to find the 
material orientation (if any) which minimizes the strain energy for a given state of stress or 
strain. New results are presented in this paper on the determination of optimal orientations 
for both general and specific types of anisotropy. 

The general problem of determining optimal orientations in anisotropic elasticity has 
been the subject of several studies in the last two decades, beginning with the work of 
Seregin and Troitskii jl6j in the context of orthotropic solids. They determined the impor- 
tant coaxiality condition: a minimum or maximum of strain energy requires that the stress 
and strain share common principal axes. The coaxiality condition was subsequently and 
independently obtained by others: first, by Rovati and Taliercio ^1] who considered 3D 
elastic materials with orthotropic and cubic symmetries (although their derivation is not 
restricted to these symmetries but is applicable to general anisotropy), and later by Cowin 
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[S]. Cowin derived the coaxiality condition independent of material symmetry considera- 
tions. He showed that the commutativity of the stress and strain is a consequence of the 
stationarity condition of the strain energy with respect to rotations of the moduh. Vianello 
[TU| provided a more formal derivation of coaxiality in linear elasticity. He used the tangent 
space of the rotation group to show that there are at least two orientations of the moduli 
that simultaneously make the energy stationary and stress and strain coaxial, a result later 
generalized to hyperelasticity jSD] (it was subsequently shown that at least three such ori- 
entations exist for both linear elasticity jT2] and for hyperelasticity ^H])- There is a slight 
distinction between the problems considered by Cowin and by Vianello that is important 
to note for our purposes 0. Thus, Cowin |S] considered stress states with fixed principal 
directions but arbitrary amplitudes, whereas Vianello assumed a specific state of stress. 
The former constraint defines a smaller set of possible elastic moduli for which coaxility 
can be attained, because it requires that optimal condition be simultaneously satisfied by 
a family of coaxial stresses. Not surprisingly, Cowin found that only materials with or- 
thotropic symmetry meet these conditions. In this paper the stress state is taken as given, 
in the same spirit as and also of Rovati and Taliercio jj^ . While the emphasis here is 
on three dimensional elasticity, the optimality problem has also been addressed within the 
context of two dimensional elasticity |12l 1131 |Hj . Cowin and Yang [J] considered a related 
but more general question of optimality with respect to Kelvin modes, rather than simply 
the freedom to orient a given material. A more extensive review of the literature can be 
found in the recent paper of Rovati and Taliercio |15| . 

It is interesting to note that the coaxiality condition has been derived in a variety of 
different ways: for particular symmetries (e.g. orthotropic) |lfij . using Lagrange multipliers 
, from general analytic considerations ^Hl , and even using the 6-dimensional eigenvector 
properties of the elasticity tensor The derivation of the coaxiality condition presented 
here differs from all these previous methods. Our starting point is a representation of the 
rotation matrix due to Mehrabadi, Cowin and Jaric |10| . This formulation also enables 
derivation of conditions for minima or maxima, in a simpler and more general form than 
that obtained by Cowin [21 . Section El begins with the problem definition and notation. 
The stationarity and min/max conditions are discussed in Sectional Specific conditions 
for particular material symmetries are derived in Section^ and we conclude in Section |S1 
by defining the strain deviation angle, a concept which could have application in practical 
circumstances in anisotropic elasticity. 

2 Problem definition and notation 

2.1 Optimal orientation of anisotropic solids 

Consider a fixed coordinate system {ei, 62, 63} coincident with the principal axes of stress. 
The stress tensor is therefore 



where cri, an and am are the principal stresses, in no particular order. Alternative expres- 
sions for the stress include the 3x3 matrix representation, 



cr = (Ti ei (g) ei + an 62 ® e2 + am 63 (g) 63, 



(1) 



CTI 
an 
am 
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and indicial^ notation, 

(T,y = (71 5ii5ji + an 6a5j2 + om S^Sj^,. (3) 

Our goal is to find the orientation or orientations which minimize the energy function for 
fixed stress 

£ = Ui-jUki Sijki ■ (4) 

Here, Sijki are the components of the fourth order comphancc tensor relative to {ei, 62, 63}. 
Think of the material as being free to be oriented in such a way that the energy £ depends 
upon the orientation of the material with respect to the fixed principal axes of the stress. 
The material moduli for stiffness and compliance are C*-"-* and 5*^°^ when aligned with the 
fixed axes. It is not necessary to specify at this stage whether or not the moduli possess any 
symmetry with respect to these axes. The main point is that the material is free to orient 
in arbitrary directions with oriented moduli C and § while the stress orientation remains 
fixed. 



2.2 Notation and tensor rotation 

Hooke's law relating stress crtj and strain Sij is 

= CijklSkl, £ij = SijklCTkl- (5) 

Here Cijki denote the components of the stiffness tensor, inverse to the compliance: CijkiSkipq = 
lijpqy where hjpq = [SipSjq + Siq6jp)/2 is the fourth order identity tensor. The rotated elas- 
ticity components could be expressed in terms of the unrotated components c^^\i and s[^],; , 
using Euler angles, for instance. The concise Voigt notation is used to represent the elements 
of the elasticity tensor in the fixed basis. Thus, the compliance is § = [Sij] , /, J = 1, 2, . . . 6, 
with S12 = S1122- <5'i6 = S1112, 5'44 S2323, etc., 
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An alternative representation for the elasticity tensor, closely related to is the 6x6 



matrix S with elements 



Sij defined as 



TST, where T 
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V21 



(7) 



^Lower case Latin suffices take on the values 1, 2, and 3, and the summation convention on repeated 
indices is assumed unless noted otherwise. 
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Explicitly, 
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This representation is useful in taking advantage of the fact that fourth order elasticity 
tensors in 3 dimensions are equivalent to second order symmetric tensor of 6 dimensions 
Pj. Similar equations follow for C and C = TCT. Define 
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then the stress strain relations lO become 



Sct. 



(10) 



Note that S and C are the matrix inverse of each other; SC = CS = I where I =diag(l, 1,1,1,1,1) 

The rotation about n, |n| = 1 by angle is defined as Q(n, 0) e 0(3), such that vectors 
(including the basis vectors) transform as v ^ v' = Qv. By considering small rotations, 
it may be readily seen that Q(n, 0) can be expressed in terms a skew symmetric tensor P 
that is linear in n. Thus, 



^(n, (/)) = P Q(n, (/)), where Pij (n) = ey^nj, , 

and hence 

Q = e*^. 

Note that Q possesses alternate well-known expressions 

Q(n,(^) = I + sin(/)P + (1 - cos0)p2 

= n ® n + cos (I — n (X) n) + sin P . 

In particular for our needs here, the small angle expansion is 

Q(n,(/)) =I + 0P + O(02). 



(11) 
(12) 

(13) 
(14) 



Under the change of basis associated with Q(n, (^), second order tensors (including stress 
and strain) transform as cr ^ cr', where 



QirQ j s ^rs ^ ij — Qijrs 



(15) 



Optimal orientation 



5 



The fourth order "rotation" tensor follows from (|15|l as 



Qijrs — 2 (^QirQjs ^" QisQjr^: 



and eqs. ifTT)) and ((TB|l imply 
where 



with Qy rs(n, 0) = /i- 



The formal solution of (|17|l . with meaning that should be clear, is 



Q 



and the small angle approximation is 



Q 



UP? ~ ^UP9 



0(02 



(16) 

(17) 
(18) 
(19) 
(20) 



Mehrabadi et al. ^U] derived an elegant expression for Q analogous to the representation 
for Q(n, 0). The key is the characteristic equation of P, {V^ + SP^ + 4P = where 
'^ijki ~ T^ijpqT^pqki etc.) which permits the exponential expression l|19() to be simplified. 
The result is most simply presented in terms of the 6x6 rotation matrix Q introduced by 
Mehrabadi et al. ^UJi and defined in the same manner as before. Thus, Q = TQT and 
P = TVT where T is defined in O and Q and V are the 6x6 Voigt matrices. Explicitly, 
P is a skew symmetric six dimensional tensor linear in n. 
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Q(n, (j)) is an orthogonal second order tensor of six dimensions, satisfying QQ-^ 
Equation ifT^ becomes 



Q(n, 



iP(n) 



Q^Q = I. 

(22) 



and has the explicit expansion |1(J| 



Q(n,< 



I+sin(7!)P + (l-cos(A)P^ 
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■sin(A(l-cos(A) fP+P'' 
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Finally, we note that fourth order tensors C transforms as C ^ C = QCQ"^. 



-p4). 
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2.3 Orientation function revisited 

Denote the matrix of rotation from the fixed to the "current" axes as Q. Thus, 

C = QC(°)Q^, S = QS(°)Q^. (24) 
Hence the objective function of Q for the stress-based energy minimization becomes 

£ = a'^Sa. (25) 



This is the starting point in the next section to derive conditions necessary for a minimum. 
It is important to emphasize the initial assumption that the stress is ahgned with the fixed 
axes, (O, or in terms of ct. 



o-i 
o-ii 

O'lII 







(26) 



This ensures that the energy varies as the material axes are rotated (if the stress were also 
rotated then the energy would be, trivially, unchanged). 



3 Stationarity and min/max conditions 
3.1 Angular derivatives of the strain energy 

Consider the energy £ of (|25() as a function of the rotation Q. A stationary value is 
obtained if £ is unchanged with respect to additional small rotations of S. This requires 
calculating the first derivative with respect to rotation angle for arbitrary rotation. The 
second derivative is needed to distinguish the stationary point as a minimum or maximum, 
or otherwise. 

The strain energy obtained by arbitrary rotation of the material about the axis n is 

f(n,(/))=^^Q(n,0)SQ^(n,(/.)5^. (27) 
The first derivative can be expressed, 

^£(n,0) = ?^[P(n)Q(n,0)SQ^(n,(/)) + Q(n,0)S^Q^(n,(/.)P^(n)]? 

= 2?^P(n)Q(n,0)SQ^(n,0)S. (28) 

This follows from eqs. H22() and (|27() . using the fact that P and Q commute. Similarly, 

the second derivative follows as 

— ^£:(n, 0) = 2^^ [P'(n)Q(n, 0)SQ^(n, <j>) + P(n)Q(n, 0)SQ^(n, (/.)P^(n)] a . (29) 
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3.2 Condition for stationary strain energy 



Assume, with no loss in generality, that the stationary orientation is at = 0. If S is at a 
stationary point, then the energy should be unchanged regardless of the axis n, or 



£ stationary <^ 



d 



= 0, V|n| = l. 



</>=o 



This becomes, using (|28)l evaluated at = 0, 

d 



dcj) 



£{n,> 



2ct P(n)S5^. 
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0=0 



We now take advantage of the fact that the stress is aligned with the fixed axes. Thus, cqs. 
(inU and ^ give 



Hence, 
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This must vanish for arbitrary direction n, hence the energy £ is stationary if 
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Let us assume for simplicity that the state of stress is triaxial, so that (Ti, (Th, a\\\ are 
distinct. The left matrix in H34|l can be removed, implying a linear condition in the stress: 



E(ai, (711, CTiii)^ = (0, 0, 0)^, 



(35) 



where E involves moduli (compliances) only 
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Thus, the energy function £ is stationary if (ctj, cth, crm) is a right null vector of the 3x3 
matrix E. Based on (IIOII and jSjl, the condition (|35|l is equivalent to the requirement that 
the off-diagonal elements of the strain vanish: 



£23 £31 — £12 — , 
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or, 
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where ei,eiieiii are the principal strains. We have therefore derived the simple but impor- 
tant general result: 

Result 1: The energy £ is stationary iff the stress and strain are coaxial. 

Equation (|37|) states that the 3-vector (cri, (Th, <^\\\) is a right null-vector of E. This 
requires as a necessary but not sufficient condition that 

detE = 0. (39) 

Consequences of this condition were explored in detail by Rovati and Taliercio |15| for 
particular material symmetries: cubic, transversely isotropic and tetragonal. A different 
approach is taken in Section^ below, where the strain energy will be minimized directly. 

While Result 1 is not new but has been derived by several authors [161 1141 |S1 1191 13 El j 
the present derivation is novel and explicit. In particular, it allows us to go further and 
obtain the condition necessary for a minimum or maximum. This is explored next. 

3.3 Condition for an energy minimum 

The second derivative of £ at the stationary point follows from (|29|l evaluated at </> = 0, 
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Each term on the right hand side will be examined in turn. Using H32|) . it follows that 
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Thus, 



4n^Fn, 



(43) 



= 



where the symmetric 3x3 matrix F has elements 
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The second derivative (|43|l must be positive for all directions n at an orientation where 
f is a local minimum. Noting that 
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i^i2 = 2545(0111 - 0]])(oi - CTiii) + ^ (o] + O]] - 2o-]n) £12, 
it follows that at a stationary point the off-diagonal elements of F become, 
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Equivalently, by pre- and post-multiplication of -^F by the diagonal matrix diag[(o-]]] 
0]])~^, (o] — oiii)^^, {ail — oi)^^], it follows that G must be positive definite, where 
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In summary, 
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This must hold for arbitrary n, |n| — 1, and again assuming that the principal stresses 
are distinct, it follows that G must be positive definite. Combined with Result 1 for the 
existence of a stationary point, this gives 

Result 2: The energy £ is a local minimum if the stress and strain are coaxial and the 
symmetric matrix G of (|48|l is positive definite. 
This can be rewritten (with obvious notation) 
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The left matrix is positive definite because of the positive definite properties of the moduli. 
The minimum condition therefore requires that this latter matrix be greater than^ the 
right hand diagonal matrix defined by the principal stresses and strains. The requirement 
that the full matrix is positive definite can be relaxed if the stationarity is restricted in 
orientation axis n. Thus, only the single scalar quantity n^Gn needs to be considered in 
the important special case of rotation about a single axis. This situation is examined in 
detail in the Appendix. 



4 Optimal orientation of particular material symme- 
tries 

4.1 Partition of the energy 

Before considering specific material symmetries, several general results are presented which 
help focus attention on the anisotropic part of the strain energy. Separate contributions to 
the energy function £ of from isotropic and anisotropic parts of the elastic moduli may 
be distinguished as follows, 

(is) I (an) ft^T\ 

= (^^(^kl S\jil + (T^jCTkl s\.^ . (51) 

The moduli are partitioned into isotropic and anisotropic parts 

t>ijkl — ^ijkl ^ ^ijkl 1 '^Vkl — Hjkl ^ '^ijkl ' W^J 

with the isotropic moduli defined by 

^ 3k' ^ 2jr ^^'^^i ^'^^ '^'^''^ ^^'^ Kijki- (53) 

^Thc matrix U is greater than the matrix V if r-^Ur > r^Vr for all nonzero r S R^. 
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Here, 
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and lijki arc the elements of the fourth order identity. The effective "bulk" and "shear" 
moduli Kg, /is and Kc, jic are obtained as 
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4/is Z 

9kc = Sc^jkiJzjki = Cii + C22 + C33 + 2Ci2 + 2Ci3 + 2C23, (57) 
3 

15/ic = T;CijkiKijki = Cii + C22 + C33 — C12 — C23 — C31 + 3C44 + 3C55 + 3C66- (58) 



Note that in general 7^ Ks and 7^ Ms- The anisotropic parts of the moduli in H52() are 
defined as the remainder after subtracting the isotropic parts, s|™; ~ Sijki — Si^^, etc. 



The energy associated with the isotropic part of the moduli becomes 
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These may be written explicitly in terms of the principal stresses, from as 
1 , 
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By definition, the scalar quantity is zero, and accordingly, the anisotropic energy 

simplifies to 

c(an) o — / (an) / / (an) /ciA\ 

It may be shown that the 3x3 matrices E and G of H36|) and (|48|l vanish for isotropic 
materials. In general, they depend upon the anisotropic part of the material moduli. 
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4.2 Materials with cubic symmetry 



In the fixed coordinate system of the principal stress axes, the elastic compliance for a 
material with cubic symmetry is 
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There are three independent moduli, k, fii and /i2, where 
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The associated fourth order tensors can be expressed succinctly using the irreducible tensor 
notation of Walpole |23 , as 
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^^jfci - = - J^Jkl, and 
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This format makes it relatively straightforward to determine the effective isotropic moduli, 
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10 



1 



-) (5J, 
M2/ 



ijkl 



2K, 



ijkl 



(70) 



The anisotropic part of the energy, (|64|l . depends only upon the dcviatoric part of the 
stress, 



£(an) 



/ / (an) 
^ij^kl Sijkl ■ 



(71) 



The reason is that the second order tensor s 



(an,0) 
ijkk 



is identically zero for cubic symmetry, and 



hence remains zero in the rotated material axes: slj'^j. = 0. The first term in H64|) therefore 
vanishes, leaving the simpler expression H71|l . The isotropic tensors Jijki and Kijki are 
unchanged under rotation, and consequently, from H59|l . (|70ll and (|71|) . 
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where Dijki is the rotated version of of^^^. In order to avoid ambiguity, let U(^j.gi^j;i denote 
the stress in the rotated coordinate system, then it foUows that 



^ij'^kl Dijkl — ""(rot)!! + '^{rot)22 + ^(rot)33 

The scalar second invariant of the dcviatoric stress is 



.'2 



J2 



(73) 



a 



<7'u + ' 



^III 



t'2 

'(rot)ll 



-'(rot)22 



'2 

''(rot)33 



'(rot)23 ^ ^"(rot)31 ^ ^'^(rot)12- C''^) 

Therefore, the function £^^^^ of H72|l is stationary when either the right member of (|73|) or 



r = cr 



'2 

(rot)23 



t'2 

'(rot)31 



are stationary. Furthermore, 

^(cx) ^ ^(cx,o) + n _ ±^ r ' 
VA^i M2/ ^ 



rot)23 



t'2 

'^(rot)12' 



'2 '2 1 

''(rot)31 + '''(T-ot)12j 1 



where £('='^'0) is the unrotated or fixed value, which follows from as 



£-(cX,0) 



1 / 1 



2VM2 Ml 



Hence, 



Ml > M2 



/^'^ I 2 / 2\ 



£;(cx) ^ £(cx,0) 



(75) 



(76) 



(77) 



(78) 



with equality when the material and stress axes are aligned. Thus, a local minimum that 
is not aligned with the stress axes occurs iff [i\ > ^2 and occurs when F of H75() achieves a 
local maximum. It will be shown below that the maximum value is ^'^ij'^ijt or equivalently, 
that is zero at the stationary point. 

'^ijkh '-'^^y part that con- 
tributes to the anisotropic strain energy is D^jj^i Dijki- Conditions for obtaining the 
stationary value of strain energy are next derived by focusing on the dependence upon 
Di-iki- The 6x6 matrix associated with the unrotated tensor D^^l^ is 



As the material axes are rotated to transform 5,-"^';'°'' 

■,(0) 



D(o) 



^3x3 03x3 
03x3 03x3 

It is convenient to split Q as follows into 3x3 matrices, 

Qi Q2 
Q3 Q4 

so that the rotated tensor D — QD("'Q^, follows from ifT^ and (|5n|l as 



(79) 



Q = 



(80) 



D 



QiQf Q1Q3 



Q3Q1 Q3Q3 



(81) 
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The term associated with the rotated energy follows from (|26|l and (|81() as 



(82) 



Thus, any stress that is a null vector of also yields the minimum or maximum value 
for ft™) of H72|) . i.e. zero. This suggests Qi as the focus of attention, and implies the 
important result: Every stress state which is a null vector of corresponds to a global 
minimum (maximum) of £ if /ii > /i2 (/i2 > yUi)- We therefore search for null vectors of 

Before deriving two alternative methods to find null vectors of in the next two 
subsections, note that the quantity l|73|) vanishes at a stationary orientation, and hence 
C(rot)ii = '^(rot)22 — ^{rot)22- Thus, the strcss iu cach of the three rotated axial directions 
is equal, a result previously obtained by Rovati and Taliercio |14lll5j . Furthermore, at the 
stationary point it may be easily shown that the following identities hold: 



Da = (7 1, Ccr = 



Ma- = 0, 



(83) 



where V^C^M are the (rotated) tensors with components Dijj^i, Lij^-i, Miji^i, respectively. 
Hence, the strain at optimal orientation is simply 



£ = (7 I 

3k 



1 



7" — a' (optimal orientation only). 
2/^1 



(84) 



This is clearly coaxial with the stress, which follows from the commutation property of 
coaxial second order tensors {ae — ecr = in this case). 

It is also worth remarking that we do not seek null vectors of the matrix E, although 
this approach is feasible and has been used to advantage by Rovati and Taliercio ^] . Some 
comments on E are in order though. The 3x3 matrix follows from eqs. (|79|l through (|81|l 
as 

E = Q3Qr, (85) 

and the condition (|39|l is satisfied if either detQi or det Q3 vanish. These can be made 
more explicit in terms of the elements of the rotation matrix. Let 



Q = 



qii 
q2i 

931 



912 913 
922 923 
932 933 



(86) 



then using the the explicit representation of the 6x6 rotation matrix from Auld |Tj or 
otherwise, the condition (|39|l implies 





9?i 


9?2 


9?3 




921931 


922932 


923933 




det Qi = 


9I1 


9I2 


9^3 


= 0, or dot Q3 = (2)3/2 


931911 


932912 


933913 


= 




9I1 


9I2 


933 




911921 


912922 


913923 





(87) 



Using the fact that the column vectors of a transformation matrix form an orthonormal 
triad, it follows that 



Q^l, 1, If = (1, 1, If, Q3 (1,1,1) 



(0, 0, 0)^ 



That is, ((Ti, (Tn, crm) = A(l, 1, 1) is a null vector of E for any A. Hence, E is not of 
full rank, implying that det E is always zero. However, this is not of interest as the null 
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space corresponds to hydrostatic stress, for which the energy function is independent of 
material orientation. The imphcations of H87|) are not considered further, and we return to 
the simpler task of finding null vectors of alone. 

Two methods for achieving the minimum energy are described, both using 

explicit forms of the rotation. The first involves a single rotation about an arbitrary axis, 
and the second is in terms of standard Euler angles. 



4.2.1 Minimum energy state v^rith a single rotation 

The range of transformations which correspond to energy minima can be obtained using 
Euler's Theorem ^ which states that any transformation matrix Q can be represented in 
the form H12|l for some axis n, |n| = 1, and angle (p. Thus, 



Q(n, 



1 — 2s'^(n2 + 2s{snin2 — cn^) 2s{snin3 + cn2) 

2s{snin2 + cna) 1 — 2s^(n| + nf) 2s{sn2n^ — cni) 
2s(snin3 — cn2) 2s{sn2n^ + crii) 1 — 2s'^{n\ + n^) 



(89) 



where c = cos^, s = sin^. and the elements of the associated Qi are determined by 
squaring each element in 189|1 . It may be shown that 



det Qi (n. 



cos 2(/)+2(2+3 cos (/>) (1— cos (fyf {n\n2+n2n\+n\n\)+&{l—cos (j))^ n'^n'^n'^ 



Note that nln^ + n^n^ 
nl = 1/3. 



^^nf < 1/3 and riln^n^ < 1/27 with equality when 



(90) 



For a given rig and angle 0, 



1 



nlnl = -{l-nl)±[-il-nl)'-g_ 



1/2 



where 



2V -o. 

- cos 2^ - 2(2 + 3 cos - cos (jjf n§(l - n§) 



2(1 - cos 0)2 [(2 + 3nl + 3 cos 0(1 - n§)] 
The null vector of Qi is such that 



(91) 
(92) 



<y'iQn + cr[iQl2 + (TniQ 



13 



0, 



(j{Q2i + <y'uQ22 + ^mQ2d. = 0' 

0. 



Using the fact that this is a dcviatoric stress, we replace crjjj 
of (HSI), to get 



o"! oil - Q33 ) + ^'\\ \Q\2-Q 



(93) 

-(jj — (tJj in the final equation 
0. 



Hence, 



- ao {q\2 - Q\-a^ ■. ^'\\ - «0 ^Q33 ~ ^31^ ' 
for arbitrary oq 7^ 0. Once again using the fact that cr(jj = — a{ — ctjj gives 



o^ni = ao Q31 ^ 



^32 



(94) 
(95) 

(96) 



Optimal orientation 



16 



In the same way, using the other equations in H93|l . three alternative expressions for the nuU 
vector are found: 





= ai ^ 


Q32 


- QI2, ' 


333 


-Q31>' 


931 


^ Q32 






Q12 


- QI2, < 


9l3 


~Qli, < 


9ii 


~ *9l2 




= asf 


Q22 


- Q22, < 


923 


-Q21, < 


921 


^ ^22 



(97) 
(98) 
(99) 



for some constants ai, a2, a^. Thus, from the first expression, with ai = 1, 
a[ = [(1 + nj)(n2 — n3)(l — cos 0) + 713 — ^2 — 4nin2n3 sin 0] (1 — cos c 



^ = {[{l-nl){nl-nl)-2nl]{l-cosc 



{[{l-nl){nj-ni) + 2nl]il 



'S0), (100) 
+ 712 - + 1 + 27ii7i27i3 sin0}(l - coscj)) - 1(101) 
+ 71^ - rig - 1 + 2nin2?i3 sin0}(l - cos 0) + 1(102) 



These equations provide us with a two parameter set of stress states, described by < 7i| < 1 
and ip. The two are independent insofar as ti^ and 77,2 of (|91|l he in (0,1). This in turn 
requires that g of 1)92(1 satisfies 



2\2 



(103) 



O<g(n^,0) < ^{l-ni) 
which defines the range of < n§ < 1 and 4>. 
4.2.2 Minimum energy using Euler angles 

The standard three Euler angles (0i, 6*2, ^3) are used to transform from {81,62,63} 
{6'i,6^,6^ = 63} {e'{ = 6i,6'2',6^j'} {e'j", 6^", 6^" = 6^'}. That is, first rotate about 
the 63 axis by 9i, then about the e'^ axis by 6*2, and finally about the 63 axis by 83. The 
transformation matrix is 

Q(^l,^2,^3) 

cos 9i cos 6*3 — sin 9i cos 62 sin 6*3 sin di cos 63 + cos 9i cos 62 sin 6*3 sin 62 sin 63 



- cos di sin 63 

sin 61 sin 62 



sin 9i cos 6*2 cos ^3 — sin 9i sin ^3 + cos 9i cos 6*2 cos ^3 sin 62 cos ^3 

— cos 01 sin 6*2 cos ^2 



and it follows from this and H87|l i that 
detQi (6*1, 6*2,613) =cos26'i cos26'2 cos26'3-i sin26'i sin26'3 cos 6*2 (3 cos 26*2 + 1) • (105) 

The condition that this vanish is equivalent to cq. (90) of Rovati and Taliercio j although 
their result is obtained in a different manner. 
Consider, for instance, 9^ = 0, for which 

detQi (6ii,6i2,0) =cos26'i cos26i2, (106) 



(104) 
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and hence there are null spaces associated with Qi(0i, 7r/4, 0) and Qi(7r/4, 6*2, 0). The null 
spaces are lines in the stress space, which follow from the simplified form of Q^, 



COS 



sin^ 9i cos^ 02 



sin^ 01 sin^ 02 



sin^ 01 cos^ 01 cos^ 02 cos^ 0i sin^ ^2 



sin^ 02 



cos 



(107) 



The possible states of deviatoric stress are: (cr{, (t{j, (t{jj) = A (0, —1, 1) if 02 = 7r/4, 6'3 = 0, 
and {a[, cr{j, crjjj) = A (cos 26*2, - cos^ 02, sin^ 02), if = 7r/4, 6*3 = 0. The first family of 
stresses correspond to a 2D elasticity problem (see Appendix): a[ = 0, crjj + ctJjj = 0, and 
it is also a null vector of Qi(0, 7r/4, ^3). The second is also a null vector of Qi(0, 6'2, 7r/4). 
Similarly, A(l, —1, 0) is a null vector of Qi(7r/4, 7r/2, ^3) and A(— cos^6'i, sin^ 6*1, cos 26*1) 
is a null vector of Qi(0i, 7r/2, 7r/4). 

Conversely, an orientation which provides a minimum in energy can be found for a given 
state of stress. Assume with no loss in generality that <t[j < < crjjj. Define the angle 02 
by 

(108) 
(109) 



tan2 02 = -^, 
^11 



then the deviatoric stress may be expressed 

(cr{, a[i, a[ii) = {a[ii - ctJi) (cos 26*2, - cos^ 02, sin^ 6*2). 



It follows from the above example that this deviatoric stress is a null vector of Qi(7r/4, 02, 0). 

It is instructive to examine further the example of eq. (|109|1 . The rotated material axes, 
denoted {e'^, e'2, 63}, are given by the columns of Q(7r/4, 02, 0): 



(110) 





1 ' 




1 " 




" " 
















— c 




c 


, e;3 = 


s 




s 




— s 




c 



where s = sin02, c = cos 6*2, or from H108() . 



'Ill 



'III 



'III 



The rotated tensor Dijki is 



V = e[i^ e\ (g) e[ (g) e'^ + (g 62 (g ® 63 + 63 (g) eg (g) 63 (g) e; 



3 ' 



and hence 



I]e;®e',(e;.^.rT'ey, 

fe=i 



(111) 

(112) 
(113) 



It may be seen by direct calculation that the three scalars e'f. ■ cr'e'f. (no sum) are identically 
zero by virtue of (|110|1 and (|lll|l . This demonstrates explicitly that 



Per' = 0, 

at the optimal orientation. The identities (|83|l follow accordingly. 



(114) 
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4.2.3 Summary for cubic symmetry 

The extreme values of the energy for cubic materials are £i and £21 where 

The fixed axes are always one of the stationary orientations, since E of H3t)|) vanishes. The 
stationary value for the unrotated axes is £2, which is the global minimum (maximum) if 
M2 > A*i (^2 < Ml)- The stationary value £1 occurs at some rotated axes, the existence 
of which is not in doubt for a material of cubic symmetry (or any material symmetry for 
that matter). The important point to note is that it is possible to explicitly determine such 
orientations. Thus, we have shown by direct construction the material orientation that 
achieves the stationary energy value £1 for any state of stress. This is a global minimum 
(maximum) if /X2 < /ii (^2 > A^i)- 

It is interesting that the expressions for the extreme values in pi5|l have the form of the 
energy for an isotropic solid, but with different shear moduli. This is evident by writing £1 
and £2 in terms of the invariants of the stress tensor: 

^ + ^f*"^' ^(tr*^)']' ^ = 1' 2- (116) 

4.2.4 Example materials 

Noting that /ii = C44 and ^2 = (cn — ci2)/2 allows us to determine the sign of (/i2 — Mi)- 
Table A.l of Musgrave provides data for c* = 2(/i2 — /ii) for a multitude of materials. 
These show c* to be negative for most elemental and engineering materials with cubic 
symmetry and different lattice structures: aluminum, nickel, copper, silver, gold (all f.c.c 
structure), iron (b.c.c), brass (f.c.c. and b.c.c), diamond, silicon, germanium (all diamond 
structure), and GaSb, InSb (both zinc-blende). Hence, for all of these cubic materials 
there exist optimal orientations of the axes that achieve the lowest energy state possible. 
Materials with positive c* include crystalline compounds of potassium, sodium and silver 
with rock-sah structure: KF, KCL, KBr, KI, NaF, NaCl, NaBr, Nal, AgCl, AgBr; plus 
caesium compounds with structure related to b.c.c. For these, the orientation associated 
with (|109|l . for instance, gives maximum strain energy. The minimum energy is achieved 
by no rotation. 

4.3 Transverse isotropy 



Materials with hexagonal symmetry, or equivalently, transverse isotropy (TI), are charac- 
terized by five moduli. In the coordinate system of the principal axes, the elements of the 
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compliance are 





11 


c.(0) 
12 


c.(0) 















c.(0) 
'-'ll 


c.(0) 
"-"13 

















c.(0) 
■^33 











§(0) ^ 






c.(0) 
J 44 










s 


Y 


M 




c(0) 
J 44 




c(0) 
'-'66 



(117) 



with S'gg^ = ^(•S^^'' — 'S'i2'')- The TI material is characterized by an axis of symmetry, 
defined by the unit vector n, which is here chosen as the axis. In general, the strain 
energy depends only upon the orientation of n with respect to the stress axes, and the 
problem is formulated as one of selecting n to minimize £. 

First note that two of the five moduli can be ascribed to the isotropic part of the 
elasticity; or conversely, an isotropic part may be subtracted from the compliance tensor 
Sijki according to (|^ . ifK^ and ifK^ . where 



:(o) 



- = 2S[°hsilh2S[',^+4S[l\ ^ = 2S[l^ + sif^sil^^2S[lh6Si"J+3Sg\ (118) 



leaving a tensor s^™; with three constants. The anisotropic compliance depends upon the 
orientation of the axis of symmetry as follows 



(an) 



aniTijUkni + h {SijUkUi + Skin^Uj) + - [SikUjUi + SiiUjUk + Sjknini + 5jin.,nk) 



-(a + 66 + 2c)Jijki - — (a + 5c)Kijki- 



(119) 



The tensors Jijki and Kijki are defined in |MJ, and the elastic constants a, b and c follow 
from l|TT7jl - fTit|) as 



a = S[f + S; 



33 



2S- 



(0) 



4.S. 



(0) 



6 = 5'- 



(0) 



(0) 
2 I 



c — 2sij 



2S, 



(0) 
66 ■ 



(120) 



The stress is, as usual, aligned with the fixed axes, so that the total strain energy follows 
from (inU), l|119|l and (jSJ, as 



where is defined in (|6()() and the extra energy term is 



(121) 



f = a (CTl?l^+CTll7T.2+(7lIl7T.3)^ + (36CT + CO'l)2CTinJ + (36CT + CCrn)2CTll7l2 + (36(T + C(TlIl)2(TlIirt3 . 

(122) 

The latter shows that the anisotropic part of the energy depends upon the TI axis 

orientation through the three parameters nf , and n|, which satisfy + n.2 + n| = 1. 
Since each of n^, and n| must be non-negative, the permissible set is the equilateral 
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triangular area A of the plane nf + ii^ + n^ — 1 ^ bounded by the hues Li : + ^ 1, 
L2 : n\ + n\ = 1 and : n\+n\ = 1. 

Consider first the possibihty that the optimal orientation lies on one of the lines i^, 
1=1,2, 3. Thus, along L3, a simple calculation using ri,| = shows that 

= a(CTii - dif {nl -Nif + 2(Tii{iha + can) - a{(7u - crif Nf on L3, (123) 

where 

a{(Tu - cTi) 

Thus, nl = Nl, is a possible optimal orientation. It must first be checked whether or not 
Nl lies in (0, 1). If this is so, and if a > 0, then an energy minimum occurs at the point 
nf = A^i, n| = 1 — iVi on L3. Similarly, a minimum occurs on Li sA = N2, n\ = \ — N2 
if N2 e (0, 1) and a > 0, and on L2 at nj = N3, n\ = 1 - iVs if iVa G (0, 1) and a > 0, where 

ao-iii + c(crii + cTiii) + 36ct 

N2 = -, r , (125) 

a(o'iii - CTiij 
,r aai + c(crin + oi) + 3&ct 

iVg = . (126) 

a[a\ - (Till) 

Now consider the possibility of the minimum occurring in the interior of A. Substitute 
n\ — \~n\ — 112 iiito (|122|l and setting the partial derivatives with respect to and rij to 
zero, yields a pair of simultaneous conditions 

(o-i - (Tin) [a [ainf + aunl + auinl) + 3bW + c{ai + am)] = 0, (127) 
(tTn - (Tin) [a [<7in\ + aun^ + cniin\) + 3ba + c{(Tu + am)] = . (128) 

These cannot be satisfied in general if the three principle stresses are distinct. We therefore 
conclude that the optimal n will lie inside A iff the stress is biaxial. For instance, if an and 
am are equal, then eq. I|127(l combined with n\ = 1 — n\ — n\ imply that E^'^^^ of H122|l 
can be expressed as a function of n\ alone, and the expression is identical in form to that 
given in H123() . Thus, the existence of a minimum inside A requires biaxiality [an = am) 
and that < iVi < 1. The associated optimal direction is not unique, but is defined by the 
cone n\ = Ni, n\+ n\ = \ — Ni (note that N^ = \ — Ni when an = am)- Again, H123|l 
indicates that the optimal orientation corresponds to a minimum (maximum) in energy if 
a > (a < 0). Thus, the sign of the elastic compliance a is crucial in determining whether 
the stationary point is a minimum or a maximum. 

These conclusions may also be confirmed by the coaxiality condition for the stress and 
strain. Thus, for arbitrary orientation 

£23 = ?^2"-3 [a{ain\ + annl + amnf) + iba + c{au + am)] , 

£31 = '^3"-! [a{<yin\ + aiinl + amnf) + "iba + c{am + <^i)] 1 (129) 

£12 = nin2 [a[ain\ + aiin\ + amn\) + iba + c[ai + an)] . 

The requirement that these simultaneously vanish is identical to the above conditions for 
the existence of the minimum inside A or along its perimeter. 

In summary, a > is a necessary condition that an energy minimum occurs at points 
inside A or along the lines Lj, j = 1, 2, 3. A minimum is achieved iff one or more of A^i, N2 
or N3 lies in (0, 1). The minimum occurs on the associated bounding line Lj or on a cone 
of directions for biaxial states of stress. Otherwise, the global energy minimum corresponds 
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to one of the vertices of A, i.e. at nf = 1 or = 1 or 713 = 1. In this default case the TI 
axis of symmetry is ahgned with one of the stress axes. These findings are in agreement 
with those of Rovati and Tahercio JHl i who stated the condition as follows: At least one of 
the principal axes of stress must lie in a plane of transverse isotropy, or alternatively, the 
TI axis must lie in a plane defined a pair of principal axes of stress. 



1 qiO) 


+ Is^ - 




c(0) 
■^13 


c(0) 1 




1/0(0) 




o(0) 
•^66 ■ 



4.4 Tetragonal symmetry 

The moduli have the same general form as in H117() , except that there is no relation between 
S'gg\ s[^i and . In this sense, tetragonal symmetry is the same as TI but with one 
additional elastic constant. The isotropic moduli are given by cq. (|118|l . and the anisotropic 
part of the compliance is 

(an) c' 
^ijki niUjUkTii + h' {SijUkTii + SkiniUj) + — {SikUjUi + SiiUjUk + Sjknini + S-jiniUk) 

1 2 

+d{piPj - qiqj){pkPi - Qkqi) - ^(a' + 66' + 2c')Jijki - ^^(a' + 5c' + 3d)Kijki- (130) 

The additional fourth order tensor as compared to TI is {piPj — qiqj){pkPi ~ QkQi), where 
{n, p, q} form an orthonormal triad. The elastic constants a', b', c' and d are 

0(0) r)MO) .MO) 
'-'66 ^'-'13 ^'-'44 1 

S[l')+Sil\ c' = 2S^^-2Sit (131) 



Compared to the TI constants a, b, c of (|12()|l . 

a' = a - d, b' = b-d, c' ^ c. (132) 

The strain energy of the tetragonal material can be split into a component similar in 
form to that for a TI material, and an additional term proportional to the constant d. The 
minimization of the TI part of the energy is as before (with a', 6', c' instead of a, 6, c), and 
depends upon the orientation of n but not p and q. The additional energy term depends 
on the deviatoric part of the stress and on these directions, 

= d[{pipj - qtqj)(J-j] ^ - ^ dcr.y-CTy-, (133) 

or, in terms of the principal stresses, 

£'^' = d {a{A[ + <j{,A', + a',A[,y ^ dn[^a[^ , (134) 

where 

K=Pl-1i (no sum). (135) 

The final term in (|134|l is independent of {n, p, q}, and it is only necessary to consider the 
quantity 

£' ^df, where / = a[A[ + crJjA^ + ct^jA^jj . (136) 
The orientation dependence is captured by the quantity /. 
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It is now demonstrated that for any given n, there is at least one set of p, q orthogonal 
to n which make £' vanish. Let p(°\q(°) be an arbitrary pair of unit vectors such that 
{n, p'*^-*, q'*^-*} form an orthonormal triad, then every possible set {n, p, q} is defined by the 
pair p, q obtained by rotation about n by angle 0: 

p(0) =cos0p("' -sin0q(°), q(0) = sin</)p(°) + cos-^qt"'. (137) 
It may then be readily verified that 

A^(0) = A^(0) cos2(/)- A^(^) sin2</-, i = l,2,3. (138) 
Equation H138|l implies that 

E'{4>) = d [/(O) COS 2(/) - /( J) sin 20] ' , (139) 

and hence 

£'{4>*) = 0, where tan 2ct>* = /(0)//( J) , (140) 

Thus, if d > 0, the situation for tetragonal symmetry is a simple addition to the TI situation: 
First find n which minimizes the TI part of the energy. Then, select the pair p, q such that 
they satisfy lfTin| . The minimum energy is then exactly that achieved by the TI part of the 
moduh (although it depends upon a', 6', c' rather than a, h and c). 

If (i < then the situation is more complicated, and the sequential minimization of first 
the TI energy and then the additional energy S' does not work, although these do define 
stationary points for the strain energy. The d— term must be taken into account when 
optimizing with respect to n, and a more complicated minimization problem is involved. 

Tetragonal symmetry represents a demarcation between the simpler higher material 
symmetries for which explicit results can be obtained, and the lower material symmetries 
which require numerical resolution, in general. Exceptions may occur, however, it is useful 
and instructive to distinguish the cubic, TI and tetragonal symmetries from those of, for 
example, monoclinic symmetry with 13 independent moduli to consider. 



5 Strain deviation angle 

5.1 Definition of the strain deviation angle 

A necessary condition for an energy minimum is that the stress and strain are coaxial. This 
is always the case in isotropic media, whereas it will be the exception rather than the rule 
under conditions of general anisotropy and arbitrary stress. According to Euler's theorem 
P] the transformation from one set of principal axes to the other can be reduced to an axis 
of rotation n, |n| = 1, and an angle of rotation (f). The stress axes have been assumed to 
coincide with the fixed axes e^ , j = 1,2,3. Let be the orthonormal axes of the strain 
tensor, then it follows that the rotation matrix is simply the matrix composed of the three 
unit vectors as columns, 

Q=[e'ie^e^]. (141) 
Let Q be represented by (|12|) . then it follows from the latter that 

eijkQjk = 2sin(/)ni. (142) 

This provides a formula to determine both the angle and the axis of rotation n. 
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The strain deviation angle is defined as the angle of rotation between the stress and 
strain axes. This angle is identically zero in isotropic materials for all stress states. In 
anisotropic materials it depends on both the material constants and the state of stress. 
However, the above analysis tells us that (j> = is a, necessary condition for energy mini- 
mization. Therefore, the magnitude of (j) provides, through a single parameter, the degree 
to which the given state of stress and material orientation are optimal. It does so without 
requiring any calculation of the energy locally or globally. It requires only that the principal 
strain axes are determined, and from those <f> can be immediately computed. 

For a given material, stress and strain, the strain deviation angle can be obtained from 
(|142|l . A more explicit method is to use the general identity for integer m: 

cosmc/) = itr(Q'") - i. (143) 

This follows from, for example, eqs. (|12() and H13() . which imply 

Q"(n, 0) = n(g) n + cosm0(I - n® n) + sinTTK^P . (144) 

For instance, m — I gives the strain deviation angle explicitly in terms of the first invariant 
of Q: 

<^ = cos-i [i(trQ-l)]. (145) 



5.2 Weak anisotropy 

Let 6^^\ j = 1, 2, 3 be the principal strains for the isotropic medium, i.e. the principal axes 

of s^^^li^ki, where aij is given by (O and/or In order to determine the strain deviation 
angle we first need to find the principal axes of strain. It is useful to express the strain as 

(146) 



(is) 



where 



7y 



(an) 



(an) (an) 
■ 422'^" + 433^111 



(147) 



It is assumed that 7 is small, so that standard perturbation analysis may be used to find 
the first correction to the directions of principal strain, {e[, 63, eg} , which to leading order 
are coincident with the stress directions: 



E 



(0) 

£ ■ - e 



(o)n-i 



7e 



■3 J "^j ' 



no sum on i. 



(148) 



„{is) 



Let E and be the isotropic Young's modulus and Poisson's ratio characterizing s-ji^^, 
(|148|l implies that, to leading order. 



then 



E 



1 + ly {aj - (Ji) 
Hence, the strain deviation angle for weak anisotropy is 

E 



1 



7?2 



I23 



7I1 



(c^I - CTii)2 ((Tii - 0-in)2 ((Tin - CTi) 



1/2 



(149) 



(150) 
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It is useful to write the stress dependence explicitly by eliminating 7, 



(an) 



(an) 
CrilS24 



(an) 
^111^34 



+ 



(JlS 



16 



ClI — (Till 

n) (an) 



(an) 



(an) 



(This, 



35 



(Tin — CTl 



(an) ^ 2-1 1/2 



(Ti - (Tn 



(151) 



This shows that the strain deviation angle depends upon the same 9 moduli that appear in 
the matrix E of eq. H36(l . 

The above formula breaks down for biaxial stress. In this case the choice of fixed axes is 
arbitrary since any orthogonal pair in the the plane spanned by the equal principal stresses 
are valid. However, the choice can be made a posteriori such that the term that would 
otherwise be singular is zero. For instance, if ai = un, then the axes ei and 62 can be 
selected such that 712 — 0. 



6 Conclusions 

The 6-dimensional notation of Mehrabadi et al. ^Uj is well suited to the problem of finding 
optimal orientations of anisotropic solids. It leads quite naturally to the main results of the 
paper, which we recapitulate: 

Result 1: The energy £ is stationary iff the stress and strain are coaxial. 

Result la: A necessary (but not sufficient) condition for this to hold is that det E = 0, 

where E is defined in H36() 

Result 2: The energy £ is a local minimum if the stress and strain are coaxial and the 
symmetric matrix G of (|48|l is positive definite. 

Result 2 provides for the first time an explicit set of conditions that must be satisfied if 
the stationary condition is a minimum or a maximum. 

Specific results are given for materials of cubic, transversely isotropic and tetragonal 
symmetries. In each case the existence of a minimum or maximum depends on the sign of 
a single elastic constant. For cubic symmetry we have several new findings. For instance, 
eqs. (|l()()(l - (|in2|l provide a two parameter set of stress states which minimize or maximize 
the strain energy if a material of cubic symmetry is rotated about an arbitrary axis n by 
angle (f) (subject to the constraint (llOSfl ). Alternatively, eqs. H108() and H109|l provide a 
means to find the optimal orientation for a given state of stress. In particular, the rotation 
of the material axes depends only upon the deviatoric stress. This demonstrates that 
the stationary (minimum or maximum) value of energy can always be achieved for cubic 
materials. Furthermore, it shows that the optimal orientation of a solid with cubic material 
symmetry is not normally aligned with the symmetry directions. 

The remainder of the new results concern the optimal orientation of TI and tetragonal 
materials, and are in general agreement with results of Rovati and Taliercio |15j obtained 
by a different procedure. However, the results obtained here are more direct and provide 
considerable insight into the nature of the optimal states for these material symmetries. 
In particular, the problem for tetragonal symmetry is very similar to that for TI, with an 
additional energy term that can be simply minimized or maximized (depending on the sign 
of the constant d of (|131|l . 

Finally, we have defined and introduced the strain deviation angle. The strain deviation 
angle is inherently anisotropic, and directly related to the problem of energy minimization 
since the angle defines the degree to which a state of stress or strain is not optimal. Future 
work will explore other consequences of this new concept. 
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Appendix: Two dimensional elasticity 

Optimal orientation for two dimensional elastic anisotropy is an important special case of 
the general problem. It was recently considered by Gea and Luo 0, and is reconsidered 
here in the context of the present theory. We will see that some of the features Gea and 
Luo obtained transfer to the 3D problem: in particular, the dependence of the minimization 
upon a single elastic constant. 

The two dimensional strain energy function is 



£{0) = Siia^ + S22(^u + 2S'i2cricrii . 



(152) 



where 5*11, S22 and S12 depend upon the angle 9 by which the material is rotated rel- 
ative to the fixed 63 axis. Consider an orthotropic material with compliance elements 
•^1? ' '^2? ' "^i^^ ' ^ee' unrotated (fixed) axes. Using the standard relations |3] for the 

transformation of the moduli, it may be shown that 



S{0)^^do {an - cri)^ ( cos 261 - A)^ 

where A is a combination of stress and moduli, 

' o-ii + (Tl \ Co 



^0, 



A = 



Co and do are moduli, 



Co — '-'11 



and bo is a constant. 



-(0) 



6o = (cTii-ai)2 4°) + -(a„ + ai)^(5l^ 



(Tll - (Tl 



c.(0) 
"^11 



S. 



(0) 

22 



25*- 



(0) 

12 



45', 



(0) 
66 ' 



(0) 



s. 



(0) 
22 



2S[l^ - f) 



(153) 



(154) 



(155) 



(156) 



It can be easily seen that the energy £ of H153|l is stationary with respect to 9 when 

cos 26* = 1, and cos 261 = -1 ^ 6* = 0, and 61 = 7r/2, (157) 

respectively. Which of these yields the smaller value for £ depends upon the sign of (IqA, or 
equivalently, the sign of (ctjj— ctj )co. Specifically, the minimum is at cos 29 =sgn[(crjj— ctj )co]. 
A third stationary value is possible if — 1 < A < 1 , and occurs at 



cos 261 = A, ^ 9 ^±9* 



(158) 



i cos ^ A. If this stationary point occurs, it follows from explicit evaluation 



where 9 

that it corresponds to a global minimum or maximum of the energy. Thus, 

£{o) ] r 

> ^bo + do (o-ii - c^i)^ X < 



£{9*) J 



sm 

cos" 9* 
0. 



(159) 



It is clear that 9 = ±9* is a repeated global minimum (maximum) if do > (do < 0). This 
is the fundamental result of Gea and Luo |H] (although their conclusion is slightly different 



Optimal orientation 



26 



since they do not start with the stress in the principal axes frame): If —1 < A < 1 and 
do > then 9 = ±0* is a repeated global minimum of £{0). Otherwise, the minimum occurs 
when cos 26* =sgn[((Tjj — o'j)co]. 

The results of Gea and Luo are now reconsidered within the context of the general 
theory applied to 2D. Based on the general theory for 3D, the 2D condition for a stationary 
orientation is 

ei2 = 0, (160) 
or, in terms of the stress, assuming for simplicity that am = 0: 

S'leCTi + S'260-ii = . (161) 

The latter is consistent with the general formulation, i.e. eq. (|34|) . under the assumption 
that (Ti and cth are distinct (if they are not distinct, then the stress-based energy function 
is constant for all material orientations). The additional condition that the stationary 
orientation is a local minimum follows from 



as 



= 4^33, (162) 

?=0 



F33 = 2{au - Tif G33 > 0. (163) 
The exact form of G33 follows from H47() with (Jm = as 

G33 = See + ^(o'li - o-i)"^ [(-Sii - 5'i2)cti + (5*21 - S'22)a'ii] . (164) 
Rearrangement gives 

G33 = J (^TII - cri)~^ [{cri + cru)c + (<ti - <Jn)d] , (165) 

where c and d are 

C=5ll-522, d = 5ii + ^22 - 2^12 - 4S'66. (166) 

The specific case of an orthotropic material is considered next. It may be shown by use 
of standard relations ^ that the combinations of moduli in H166() transform according to 

c{e) = CO cos 261, d{9) = do cos 46*, (167) 

where Cq = c(0) and do = d{Q) are the same moduli defined in H155|l . Also, 

5'i6 = -i(co + docos26') sin26', (168) 

526 = -i(co - docos26') sin26'. (169) 

Hence, 

£12 = -i sin 261 [(cti + crii)co + {cri - crii)do cos 26*] . (170) 
The strain £12 vanishes if 

sin26'==0, or cos20 = A. (171) 

Thus, the stationary points are 9 — 0, 7r/2 and ±9* where cos 20* — A, in agreement with 
the results of Gea and Luo |H| ■ 
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Using the same notation, (|165|l becomes 




(172) 



In particular, if — 1 < A < 1, then 



G33{0*) = 2^0 sin2 29*. 



(173) 



This imphes that = ±0* is a local minimum of £ iff do is positive. The general analysis 
for 3D optimal orientation does not provide an explicit statement about global minima. In 
order to show that it is a global maximum one must compare the value of the energy £ at 
= ±9* with its value at the other local minimum, as done in (|159|l . 
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